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摘要 :不 同 影响 因子 对 土壤 有 机 质 含 量 的 影响 存在 尺度 依赖 性 。 以 太原 盆地 土壤 有 机 质 为 研究 对 象 , 寺 盆地 上 中、 下 部 分 别 设 
置 采样 带 ,应 用 多 元 经 验 模 态 分 解 (im multivariate empirical mode decomposition ，MEMD ) ,分 析 了 盆地 内 不 同 部 位 土壤 有 机 质 与 影 
向 因子 (高程 .坡度 ,地形 湿度 指数 .土壤 容重 砂粒. 坏 粒 、 黏 粒 和 光谱 主 份 等 ) 在 表征 民 麻 的 相关 性 ,并 预测 了 采样 尺度 上 土壤 
有 机 质 含量 , 旨 在 研究 黄土 高 原 倪 地 区 内 土壤 有 机 质 与 相关 因子 的 空间 多 尺度 美 系 。 研 究 结果 表明 :(1) 利 用 MEMD 法 可 将 分 
地 内 不 同 部 位 处 的 土壤 有 机 质 空 间 序列 分 解 为 不 同 表 征 尺 度 , 例 地 上 、 中 和 下 部 的 表征 尺度 分 别 为 6.8 和 7 个。 研究 区 域内 ， 
尺度 约 1000 m 处 是 土壤 有 机 质 的 主要 表征 尺度 , 且 盆 地 内 垂直 河流 方向 的 有 栅 质 序列 主要 表征 尺度 沿 河流 方向 表现 分 散 。 
(2) 土壤 有 机 质 和 影响 因子 的 空间 多 尺度 关系 表明 ,高程 与 土壤 有 机 质 的 关系 主要 表现 在 大 尺度 ,而 坡度 .地形 湿 度 指 数 与 盆 
地 中 .下 部 土壤 有 机 质 的 关系 较 明 显 。 土 壤 容重 与 有 机 质 在 不 同位 置 的 不 同 表 征 尺 度 存在 显著 差异 。 土 壤 质 地 中 , 壤 粒 含量 与 
有 机 质 的 多 尺度 关系 最 为 明显 。 光 谱 主 份 1 在 全 部 样 带 中 所 有 睹 征 尺 虚 下 均 与 有 机 质 显著 相关 。(3) 采 用 MEMD 法 对 有 机 质 
的 预测 精度 高 于 基于 原始 数据 的 逐步 多 元 回归 结果 。 综 上 ,研究 结果 阿 为 黄土 高 原 盆 地 区 内 土壤 数字 制图 .土壤 田 块 的 合理 设 
计 与 有 机 质 的 精确 预测 提供 理论 依据 。 

关键 词 :多 元 经 验 模 态 分 解 ; 本 征 模 函 数 ;多 尺度 ;土壤 有 机 质 ; 影 响 因 子 
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Abstract: Soil organic matter\( SOM ) content is affected by a variety of factors that operate at different scales. The purpose 
of the present study was to investigate the multi-scale spatial relationships between SOM content and influencing factors, 
including elév¥ation, slope, topographic wetness index, soil bulk density, sand content, silt content, clay content, and soil 
spectral componerits. Soil samples were collected from the Taiyuan basin, a typical area of the Chinese Loess Plateau, by 
establishing sampling transects at the upper, middle, and lower parts of the basin. At different locations in the basin, the 
multisscale correlations of SOM content with the influencing factors were analyzed, and the SOM contents at the sampling 
scale Were predicted using the multivariate empirical mode decomposition. The results showed that (1) the multivariate 
empirical mode decomposition method could separate the transects of SOM content into six, eight, and seven intrinsic mode 


functions for the upper, middle, and lower parts of the basin, respectively. A scale of 1000 m was the main representative 
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scale for SOM content in the entire basin, and the number of representative scales increased along the river. (2 ) In contrast 
to the correlation at the sampling scale, the multi-scale spatial correlation between SOM content and the influencing factors 
revealed that the correlation of elevation and SOM content was dominant at larger scales. Meanwhile, the correlation of slope 
and SOM content was significant for the middle part of the basin, and the relationship between the topographic wetness 
index and SOM content was significant for both the middle and lower parts. However, the correlation between soil bulk 
density and SOM content was much more complex and differed at the various scales and locations; in the upper part of,the 
basin, the relationship between the silt and SOM contents was more apparent than that of the sand and clay contents. In 
addition ，spectral component 1 was significantly correlated with SOM content in the entire basin. (3) The multivariate 
empirical mode _ decomposition method was more accurate at predicting SOM content than the stepwise_ multiple linear 
regression. Therefore, taken together, the results of the present study provide a basis for soil digital mappin®, dimension 


design of farmland, and SOM content prediction on the Chinese Loess Plateau. 


Key Words: multivariate empirical mode decomposition; intrinsic mode function; multi-scalel soil organic matter; 


influencing factors 


土壤 有 机 质 含量 不 仅 是 肥力 的 重要 指标 '! ,也 是 全 球 碳 循环 中 重要 的 源 和 汇 '*“”。 前 人 已 对 不 同 尺度 
土壤 有 机 质 的 空间 变异 性 进行 了 大 量 人 研究 ,证 明了 在 特定 尺度 下 ,者 二 的 异 质 性 与 外 部 环境 因素 共 WD 
s 间 分 布 "| 。 因 此 ,建立 土壤 有 机 质 与 影响 因子 的 准确 美 系 成 为 间接 理解 有 机 质 空间 变异 性 的 方法 

。 然 而 ,由 于 不 同 尺 度 下 的 土壤 有 机 质 含 量 与 影响 因 予 的 空间 关系 具有 和 较 大 差异 性 ”9 , 当 某 一 因子 在 

a 影响 土壤 有 机 质 空 间 分 布 时 ,可 能 造成 该 因 守 与 中 经 有 机 质 的 非 线 性 关系 。 因 此 ,应 用 传统 统 
计 方 法 建立 土壤 有 机 质 预测 模型 时 ,模型 的 精确 度 与 可 靠 性 可 能 会 受到 较 大 影响 。 

当前 ,在 环境 因子 与 土壤 属性 的 空间 多 尺度 关系 的 研究 中 ,多 采用 Pearson 相关 分 析 '"| 小 波 分 析 "* 沾 和 
多 维 分 析 ': 9 等 方法 。 应 用 该 类 方法 的 前 提 是 假 渴 相关 环境 因子 对 土壤 属性 空间 分 布 的 影响 为 线性 且 平 稳 ， 
但 此 类 假设 可 能 与 土壤 属性 的 实际 空间 分 布 不 相符 。 多 元 经 验 模 态 分 解 (multivariate empirical mode 
decomposition ，MEMD ) 由 Huang 等 首 软 提出 ,是 与 小 波 分 析 相 对 应 的 另 一 种 多 尺度 分 析 方 法 ,可 根据 数据 本 
身 的 尺度 特征 经 验 性 的 将 空间 数据 分 解 在 多 个 表征 尺度 上 ,避免 了 线性 和 平稳 性 的 假设 。 尽 管 MEMD 法 
具有 明显 优势 ,但 尚未 广泛 应 用 于 纪 壤 属性 尺度 问题 的 相关 研究 中 。 

黄土 高 原由 第 四 纪 冰 袁 厚 层 黄 半 堆积 而 成 ,是 全 球 水 土 流失 最 严重 的 地 区 之 一 '2 。 黄 土 高 原 盆 地 内 多 
平原 与 丘陵 ,海拔 高度 相对 较 低 , 温 湿 环 境 适宜 ,有 利于 农业 和 生产. 。 由 于 土地 利用 形式 存在 多 样 性 及 植被 
覆盖 的 斑 块 化 与 破碎 化 ,造成 了 本 区 域内 土壤 有 机 质 分 布 的 无 序 性 及 影响 因子 对 有 机 质 作用 的 非 线 性 "1 。 
因此 ,研究 该 区 域 土壤 有 机 质 与 相关 因子 间 的 空间 多 尺度 关系 可 为 有 机 质 管 理 提供 理论 依据 。 鉴 于 MEMD 
法 可 用 于 空间 数据 的 非 线性 和 非 平 稳 的 空间 多 尺度 分 析 , 因 此 本 研究 假设 该 方法 也 可 用 于 土壤 有 机 质 与 相关 
因 拖 间 的 空间 多 尺度 研究 中 。 本 文 以 山西 省 太原 盆地 为 研究 区 域 ,分 析 盆 地 内 不 同 部 位 土壤 有 机 质 的 表征 尺 
度 ,以 及 土壤 有 机 质 与 影响 因子 的 空间 多 尺度 关系 ,并 实现 有 机 质 含量 在 采样 尺度 上 的 预测 。 


了 0 研究 材料 与 方法 


1.1 研究 区 概况 

太原 命 地 地 处 黄土 高 原 中 东部 .山西 省 中 部 (36°55' 一 38°12'N,111°42' 一 113°02'E) , 属 典型 的 黄土 集中 
分 布地 带 。 盆 地 南北 长 约 150 km ,东西 宽 为 12 一 60 km, 区 域 包 括 汾 河流 域 中 游 ,总 面积 达 6159 km2 ,人口 数 
量 占 全 省 近 1/3。 该 区 气候 属 暖 温带 大 陆 性 季风 气候 类 型 ,年 日 照 总 时 数 为 2360 一 2796 bh, 年 平均 降水 量 为 
420 一 457 mm ,总 降水 量 由 南 到 北 逐 渐 减少 。 境 内 平均 海拔 800 一 900 m, 东 、 西 .北三 面 环 山 , 中 部 为 冲积 平 
原 ,年 平均 气温 7.8 一 10.3% 。 中 部 平原 以 冲 洪 积 亚 砂 粘土 质 黄土 为 主 , 边 山 丘 陵 以 风 积 厚 层 亚 砂 黄 十 为 主 。 
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境内 土壤 类 型 主要 有 潮 土 盐 化 潮 土 .石灰 性 褐 土 和 褐 土 性 土 等 四 类 ,主要 农作物 为 冬小麦 和 玉米 ,是 山西 省 
主要 农业 生产 区 。 
1.2 样品 的 采集 与 处 理 

根据 海拔 与 汾 河流 向 , 沿 北 东 -南西 方向 将 盆地 划分 为 上 .中 .下 三 部 位 。 结 合 野 外 调查 和 相关 图 件 分 析 ， 
避 开 大 型 建设 用 地 后 ,在 盆地 不 同位 置 设置 3 条 采样 带 ( 带 1 海拔 :742 一 894 m; 带 2 海拔 :723 一 807 m; 带 3 
海拔 :707 一 1002 m) ,每 条 样 带 长 约 42 km, 以 330 m 间隔 设置 采样 点 。 若 某 一 样 点 落 于 建设 用 地 或 道路 等 非 
耕地 上 , 则 在 离 该 点 最 近 的 耕地 上 取样 ,并 尽量 使 所 有 样 点 在 一 条 直线 上 。 采 样 带 1.2 和 3 分 别 包含 2 128 
和 134 个 样 点 ,共计 383 个 (图 1)。 
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1 研究 区 地 理 位 置 及 采样 点 分 布 
Fig.l Geographical location ostudy area and sample points distribution 


土 样 采 于 2016 年 3 月 22 一 31 日 。 采 样 时 首先 利用 GPS 查找 样 点 位 置 ,并 采用 “S” 形 布点 法 进行 样品 
采集 。 使 用 螺旋 取 土 钻 钻 取 0 一 20 ce 后 耕 层 土壤 样品 5 个 ,混合 后 作为 本 采样 点 样品 。 样 品 混合 均匀 后 经 风 
干 . 磨 碎 . 过 2 mm 孔 筛 后 分 为 两 份 , 汾 询 用 于 土壤 光谱 和 土壤 有 机 质 、 质 地 测试 。 使 用 环 刀 (高 5 ecm, 体积 
100 cm ) 在 “S" 形 样 点 中 心 处 采集 表层 原状 土 样 , 烘 王后 (105% ,10 hb) 测定 土壤 容重 。 采 用 重 铬 酸 钾 -外 加 热 
法 测定 土壤 有 机 质 含量 :' 光 洒 用 沉降 法 测定 土壤 质地 '' ,采用 地 物 光 谱 仪 (美国 ASD FieldSpec3 ) 测 定 土壤 
光谱 。 其 中 ,光谱 范围 为 850 一 2500 nm ,数据 重 采样 间隔 为 1 nmt"]。 

1.3 MEMD 原理 

MEMD 是 经 验 模 态 分 解 (empirical mode decomposition，EMD ) 的 扩展 ,可 直接 将 获取 的 数据 应 用 于 空间 
域 , 并 根据 室 间 数据 特征 ,将 其 分 解 为 多 个 空间 序列 。MEMD 对 数据 的 分 解 基于 如 下 假设 , 即 在 给 定 的 空间 
内 ,空间 数据 存在 简单 的 不 同 频率 相互 释 加 的 震荡 模式 。 自 然 界 中 ,事物 整体 变异 性 受 多 种 过 程 影响 ,并 
在 不 同 尺度 以 不 同 强度 发 生 避 。 其 中 ,同一 尺度 发 生 的 过 程 可 以 分 解 为 相同 的 本 征 模 函数 (Intrinsic Mode 
Function, TIMF)。 定 义 IMF 需 满足 以 下 条 件 :1) 可 以 是 线性 或 非 线 性 ,该 IMF 在 整个 数据 范围 内 ,局 部 极 值 点 
和 过 零点 的 个 数 必 须 相 等 或 最 多 相差 一 个 。 其 中 ,局 部 极 值 点 表示 一 个 函数 可 以 提取 的 最 小 值 或 最 大 值 ,过 
零点 表示 函数 改变 其 符号 的 点 。2) 该 数据 的 震荡 对 局 部 平均 值 对 称 , 即 在 任 一 数据 点 处 由 局 部 极 大 值 和 极 
小 值 定义 的 平均 包 络 值 为 0。 上 包 络 线 和 下 包 络 线 由 通过 分 别 样 条 插值 全 部 局 部 最 大 值 或 最 小 值 产生 。 根 
据 该 定义 ,IMF 可 通过 “筛选 (sifting)” 过 程 分 解 任意 函数 后 获取 ,IMF 的 方差 贡献 百分比 (%) 可 通过 下 列 公 


个 下 
方差 贡献 率 (%) = 竺 个 IME 的 广大。 x 100 (1) 
> (IMF 方差 ) + 残 差 方差 
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单个 IMF 对 整体 方差 的 贡献 率 决定 了 每 个 IMF 的 相对 重要 性 ,因此 ,也 可 决定 在 不 同 尺 度 上 发 生 过 程 的 
相对 重要 性 '” , 即 方差 分 解法 。 通 过 计算 每 个 IMF 震荡 的 平均 次 数 ,可 确定 对 应 IMF 的 平均 尺度 。 然 而 ,由 
于 土壤 属性 与 环境 因子 之 间 存 在 非 线性 关系 ,使 得 局 部 范围 内 无 法 预测 其 周期 或 尺度 。 瞬 时 尺度 可 提供 每 个 
IMF 所 代表 的 尺度 范围 ,并 通过 Hilbert 变换 得 到 的 瞬时 频率 来 获 了 到。EMD 可 将 空间 数据 分 解 为 不 同 的 表征 
尺度 ,而 Hilbert 谱 分 析 可 计算 它们 在 每 个 尺度 和 位 置 的 能 量 与 频率 ,便于 从 不 同 尺度 的 数据 中 提取 不 同位 置 
的 频率 。 

应 用 MEMD 时 要 求 计算 空间 数据 序列 的 局 部 平均 值 ,而 对 于 多 个 空间 数据 序列 可 能 无 法 直接 定 戏 最 大 
值 和 最 小 值 。 为 解决 该 问题 ,Rehman 和 Mandicl21 建 议 ,在 N 维 空间 内 通过 不 同方 向 的 投影 创建 多 个 IN_ 维 包 
络 , 取 其 均值 可 作为 局 部 平均 值 。 

假设 代表 土壤 有 机 质 与 影响 因子 的 多 维和 撩 量 数据 VGs)= jw(s) ,os)，…， ww(s)| 是 空间 数据 s 的 函 
数 , XX% = | xt ,x ,…,x*| 是 角度 矢量 = (0:,0:,…,0*_ 1) (k=1,2,…, ,KK 是 所 有 力 向 的 个 数 ) 定 义 的 沿 
不 同方 向 的 矢量 数据 , 非 线 性 的 N 个 空间 数据 的 IMF 可 通过 MEMD 的 下 列 算 法 实现 2) 。 

(1) 产 生 一 个 合适 的 方向 矢量 X 数据 集 ; 

(2) 计算 空间 序列 V(s) 沿 给 定 方向 X% 的 投影 p*(s); 

(3) 识 别 投影 p*(s) 最 大 值 出 现 的 位 置 s%; 

(4) 通 过 搬 值 [s… ,TY(s%) ] 获取 多 元 包 络 曲线 e*(s) ; 

(5) 通 过 下 式 计 算 包 络 曲 线 的 均值 W(s) ; 


M(s) -这 ew(s) (2) 


(6) 利 用 D(s)= V(s)-M(s) 提 取 D(s)。 寿 D(s) 满 是 终止 迭代 准则 , 则 从 步骤 (1) 开 始 重 复 以 上 步 又 ,并 
计算 残 差 Y(s)-D(s)。 和 否则 ,从 步骤 (2) 开 始 重复 计算 D(s)。 残 差 变 为 单调 函数 且 再 无 法 提取 IMF 时 ,分 解 
过 程 终 止 。 该 残 差 可 表明 原始 数据 的 变化 趋势 
1.4 数据 预 处 理 

去 除 土壤 光谱 数据 中 噪声 影响 较 拓 的 850 一 399 和 2451 一 2500 nm 边缘 波段 ,利用 MATLAB 将 可 见 光 - 近 
红外 波段 (401 一 2450 nm) 2050 个 竺 壤 这 详 数 据 进行 主 成 分 压缩 ,并 将 压缩 后 前 两 个 主 份 ( 占 整 体 方差 的 95% 
以 上 ) 选 作 综 合影 响 因子 。 

将 土壤 有 机 质 与 地 形 因 让 (高程 .坡度 和 地 形 湿度 指数 ) 物理 性 状 ( 容重、 砂粒、 培 粒 和 黏 粒 含量 ) .光谱 
主 份 数据 等 组 成 多 元 数据 序列 ,利用 MEMD 法 将 各 样 带 数据 序列 分 解 为 不 同 IMF。 根 据 MEMD 算法 ,通过 多 
元 序列 数据 中 的 N 个 YMF 工 同 震荡 模式 来 识别 每 个 “共同 尺度 "5 。 在 不 同 的 数据 源 中 ,MEMD 将 相似 的 尺 
度 归 为 一 组 ,代表 相关 发 生 过 程 的 实际 尺度 。 

计算 答 样 带 在 采样 尺度 与 单个 IMF( 表 征 尺 度 ) 和 残 差 下 的 土壤 有 机 质 与 地 形 因 子 .土壤 物理 性 状 、 土 壤 
光谱 主 份 等 因 淘 的 Pearson 相关 系数 。 采 用 逐步 多 元 回归 模型 ,利用 地 形 .土壤 物理 性 状 和 光谱 主 份 等 被 分 
解 的 争 IMF 和 残 差 ,预测 土壤 有 机 质 在 对 应 尺度 和 残 差 处 的 含量 值 。 最 后 ,通过 单个 IMF 和 残 差 预测 出 的 对 
应 土壤 有 机 质 值 ,利用 如 下 公式 估 测 采样 尺度 上 有 机 质 含 量 : 


SN = > IMF’? + R? (3) 
其 中 SN" 为 有 机 质 在 采样 尺度 上 的 预测 值 ,m 为 IMF 个 数 , IMF? 为 有 机 质 在 每 个 IMF 上 的 预测 值 , R* 为 
有 机 质 残 差 的 预测 值 。 
采用 由 实测 值 和 预测 值 计 算出 的 决定 系数 (coefficient of determination， 尼 ) 、 均 方 根 误差 (root mean 
squared error，RMSE ) .标准化 均 方 根 误差 (normalized root mean square error，NRMSE ) 和 相对 预测 偏差 
(residual prediction deviation, RPD ) 评价 土壤 有 机 质 预测 精度 ,各 参数 计算 公式 如 下 : 
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1 
( yeasd = i ) 2 
i=1 


R =1 -一 一 (4) 
> ( YMeasured a 也 ) 和 
i=1 


1 n 2 
RMSE = 1 CY 六 eteds) (5) 
nN i=1 
RMSE 
NRMSE = 一 一 (6 
SD 
PD = -一 一 7 
RPD = RMSE 0 
其 中 为 土壤 有 机 质 采样 点 数 ,Y™ 和 YY 分 别 为 采样 点 i 处 土壤 有 机 质 的 实测 慎 和 预测 值 , 7 为 


土壤 有 机 质 实测 值 的 平均 值 ,5D 为 土壤 有 机 质 实 测 值 的 标准 差 。 小 波 能 量 谱 代 表土 壤 有 机 质 在 空间 (位 置 ) 
和 频率 (尺度 ) 处 的 方差 ,也 可 用 于 评价 不 同方 法 对 土壤 有 机 质 的 预测 能 力 。 本 研究 中 相关 处 理 在 MATLAB 
二 中 编程 实现 。 


QO。 2。 结果 与 分 析 


2.1 采样 尺度 上 土壤 有 机 质 与 影响 因子 的 关系 

©O 表 1 为 利用 Pearson's 相关 分 析 法 对 土壤 有 机 质 与 地 形 因子 , 填 壤 物理 性 状 、 光 谱 主 份 等 影响 因子 的 线性 
本 ”相关 度 进行 分 析 的 结果 。 在 样 带 1 上 ,有 机 质 含量 与 高 程 .容重 和 郊游 主 份 1 显著 相关 ;在 样 带 2 上 ,与 高 程 、 
5 ”坡度 ,容重 ,砂粒 \ 壤 粒 、 笑 粒 和 光谱 主 份 1 显著 相关 ;在 样 带 3 和 EE 高程. 容重、 砂粒、 壤 粒 和 光谱 主 份 1 显著 
一 ”相关 。 在 全 部 样 带 上 ,有 机 质 与 高 程 . 容 重 、 砂 粒 壤 粒 和 先 谱 主 份 1 显著 相关 ,表明 在 采样 尺度 上 研究 区 内 土 
之 壤 有 机 质 主要 受 这 五 个 影响 因子 的 作用 ,相关 局 部 尺度 的 发 生 过 程 (如 生物 活动 .耕作 措施 等 ) 扰乱 了 其 他 影 
二 ” ， 响 因 子 对 土壤 有 机 质 的 作用 ,因此 表现 为 无 显著 相关 性 。 土 壤 有 机 质 与 物理 性 状 的 相关 性 顺序 为 样 带 2> 样 


一 带 3> 样 带 1。 样 带 2 中 ,地 形 因子 ,容重 . 顺 地 对 斗 壤 有 机 质 含量 的 影响 最 强 。 这 些 结果 表明 ,在 采样 太 度 上 ， 
“已 地 形 因 子 与 土壤 物理 性 状 对 盆地 中 部 条 培 有 机 质 含量 的 影响 最 强 。 


表 1 梁 样 尺度 上 土壤 有 机 质 与 影响 因子 的 相关 系数 


gp Table 1 Correlation coefficients between SOM and influencing factors at the sampling scale 
一 和 湿度 指 妆 容重 、 ed 人 ee 
样 带 高 多 放生 粒 光谱 主 份 1 。 光谱 主 份 2 
Transect Elevation Slope . Mi y Sand Silt Clay PC1 PC2 
index density 
样 带 1 Transect 1 -023 —0.01 0.11 一 0.22“ —0.04 0.11 一 0.07 0.67 一 一 0.08 
样 带 2 Transect 2 一 0.25 0.22* 0.03 -0.29** -0.52** 0.51”* 0.27”* 0.62”** -0.13 
样 带 3 Transect 3 一 0.21” -0.12 一 0.04 一 0.25 一 0.17 0.26” -0.07 0.67 一 一 0.09 
全 部 样 点 区 ll sample8 一 0.22 -0.02 0.06 一 0.22 -0.26” 0.32” 0.05 0.61 一 0.08 


吕 显著 性 来 平 为 PR0:05; * * 显著 性 水 平 为 P <0.01 


如 表 2 所 示 ,本 研究 采用 逐步 多 元 回归 法 构建 了 土壤 有 机 质 预测 模型 。 各 相关 因子 ( 样 带 1 中 为 壤 粒 、 光 
谱 主 份 1 和 2 , 样 带 2 中 为 坡度 .容重 . 壤 粒 和 光谱 主 份 1, 样 带 3 中 为 高 程 .容重 砂粒、 黏 粒 、 光 谱 主 份 1 和 2) 
对 样 带 1.2 和 3 中 土壤 有 机 质 含量 的 预测 能 力 分 别 达到 了 52% 56% 和 54% 。 
2.2 土壤 有 机 质 与 影响 因子 的 空间 多 尺度 关系 

利用 MEMD 法 将 土壤 有 机 质 与 影响 因子 的 多 元 数据 序列 分 解 为 不 同 的 IMF。 土 壤 有 机 质 的 IMF 如 图 2 
所 示 。 随 着 IMF 增 大 ,土壤 有 机 质 空间 序列 的 震荡 周期 变 长 ,表明 土壤 有 机 质 的 表征 尺度 随 IMF 变 大 而 增 
加 。 同 时 ,各 样 带 在 IMF1 处 的 震荡 幅度 均 较 大 ,表明 在 该 尺度 处 的 空间 变异 性 较 大 。 土 壤 有 机 质 的 残 差 序 
列表 明了 原始 数据 的 变换 趋势 , 即 样 带 1 中 土壤 有 机 质 的 变化 趋势 比较 明显 , 且 与 样 带 2 一 起 呈 逐 渐 减 小 的 
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趋势 ,而 样 带 3 呈 先 稳定 后 增 大 的 趋势 。 各 样 华 上 ,土壤 有 机 质 与 影响 因子 均 具 有 相同 的 IMF 个 数 与 类 似 的 
震荡 幅度 。 因 此 ,本 研究 利用 Hilbert 变换 识别 各 样 带 不 同 因 子 的 IMF 对 应 的 空间 尺度 ,并 将 所 有 因子 的 尺度 
取 其 均值 ,获取 该 样 带 上 每 个 IMF 的 对 应 尺度 ( 表 3)。 样 带 中 各 因子 的 IMF 对 应 尺度 的 变异 系数 随 IMF 变 
大 而 增 大 ,表明 在 小 尺度 处 土壤 有 机 质 和 影响 因子 具有 相近 的 表征 尺度 。 样 带 2 和 3 中 ,1 一 4 IMF 尺度 较 接 
近 , 即 在 < 5000 m 尺度 的 发 生 过 程 约 在 表征 尺度 960 .1500 .2600 和 4400 m 处 ;5 一 8 IMF 尺度 差异 较 大 , 即 样 
带 2 发 生 过 程 的 表征 尺度 约 为 8573 10901 .11591 和 23659 m, 而 样 带 3 约 为 6753、11806 和 12292 m。 该 结果 
表明 ,盆地 中 、 下 部 在 小 尺度 处 的 表征 尺度 较 接 近 ; 随 着 尺度 的 变 大 ,表征 尺度 的 差异 性 增 大 。 另 外 , 盆地 让 部 
的 表征 尺度 与 贫 地 中 部 、 下 部 的 差异 较 大 。 


表 2 基于 影响 因子 的 土壤 有 机 质 的 逐步 多 元 回归 预测 模型 


Table 2 The multiple stepwise linear regression models of SOM with the influencing factors 


样 带 Transect 公式 Function 下 R? 
样 带 1 Transect 1 3.70+16.55( 0.27)Silt+4.22(0.72) PC1-7.76( 0.16) PC2 427 0.52 
样 带 2 Transect 2 2.54+0.44( 0.17) Slope-8.94( 0.20) BD+13.79( 0.30)Silt+3.12(0.50) PC1 38.4 0.56 
样 带 3 Transect 3 49.06—0.02( 0.13) Elevation—11.31(0.16) BD-17.26(0.26) Sand-1.15(0.21) Clay+ 


3.57(0.60) PC1-5.41(0.11)PC2 2 Ds 
Elevation: 高程 ;Slope: 坡度 ;TWI: 地 形 湿度 指数 ， topographic wetness index;BD : 容重 ,Bulk density ; Sand: WSI. 壤 粒 ;Clay: 黏 粒 ;PC1 :光谱 
主 份 1;PC2: 光 谱 主 份 2; 括 号 中 的 数据 表示 标准 偏 回归 系数 


IMF1 


IMF2 


IMF3 


IMF4 


IMFS 


IMF6 


本 征 模 函数 Intrinsic mode function 


IMF7 
心 


残 差 IMF8 
写 


0 20 40 60 80 100 120 
采样 点 位 置 Sampling locations 


2 三 条 样 带 中 土壤 有 机 质 被 分 解 的 IMF 和 残 差 
Fig.2 Separated IMFs and residues for SOM at the three transects 


样 带 中 有 机 质 与 影响 因子 的 IMF 方差 百分比 如 表 4 所 示 。 样 带 1 中 ,土壤 有 机 质 方差 百分比 主要 分 布 在 
IMF1 和 2 处 (18.62% 和 19.41% ) ,表明 盆地 上 部 的 土壤 有 机 质变 异性 主要 在 1011 和 1725 m 尺度 处 。 样 带 2 
中 ,其 方差 百分比 主要 分 布 在 IMF1 和 5 处 (30.61% 和 16.91% ) ,表明 盆地 中 部 的 土壤 有 机 质变 异性 主要 在 
982 和 8573 m 尺度 处 。 样 带 3 中 ,其 方差 百分比 主要 分 布 在 IMF1.5 和 6 处 (17.62% 12.53% 和 20.46% ) , 表 
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明 该 处 土壤 有 机 质变 异性 主要 表现 在 960 .6753 和 11806 m 尺度 处 。 总 之 ,整个 区 域内 尺度 约 1000 m 均 表 现 
F 尺 度 是 小 尺度 , 即 < 2000 m; 盆 地 中 部 
的 是 小 尺度 和 中 尺度 , 即 1000、8000 m; 盆 地 下 部 的 是 小 尺度 .中 尺度 和 大 尺度 , 即 1000 .7000 .12000 m。 这 些 


为 土壤 有 机 质 的 主要 表征 尺度 。 同 时 ,盆地 上 部 土壤 有 机 质 的 主要 表 行 


结果 表明 ,该 区 域内 垂直 于 汾 河流 向 的 土壤 有 机 质 主要 表 行 


E 尺 度 沿 河流 方向 表现 分 散 。 


表 3 三 个 样 带 中 土壤 有 机 质 IMF 的 对 应 尺度 /m 
Table 3 Scales of IMFs for SOM at the three transects 
样 带 Transect IMF!1 IMF2 IMF3 IMF4 IMFS IMF6 IMF7 JWIFS 
样 带 1 Transect 1 1011(2) 1725(3) 3078(3) 5406(7) 9535(10) ”12375 (14) 
样 带 2 Transect 2 982(5) 1530( 4) 2588(6) 4487(10) 8573(14) 10901 (17) 11591 (19) 23659(3) 
样 带 3 Transect 3 960(3) 1504(3) 2604(5) 4429(9) 6753(12) 11806 (17) 12292 (她 ) 
括号 中 的 数据 表示 土壤 有 机 质 和 影响 因子 IMF 对 应 尺度 的 变异 系数 值 (% ) 
表 4 三 条 样 带 中 土壤 有 机 质 的 每 个 IMF 和 残 差 占 原 始 数据 方差 的 百分比 
Table 4 Percentage of variance explained by each IMF and residue for SOM at the three transects 
a -0g IMF1 IMF2 IMF3 IMF4 IMFS IMEF6 IMF7 IMF8 Ra 
样 带 1 有 机 质 18.62 19.41 5.12 6.56 2. 驳 8.50 一 一 24.19 
Transect 1 高 程 0.26 1.37 1.69 3.94 2.67 2.61 一 一 75.93 
坡度 39.04 18.52 10.59 13.34 4.31 3.07 一 一 6.08 
湿度 指数 25.12 15.00 7.59 10.60 5.34 8.62 一 一 7.18 
本 容重 25.12 15.00 7.59 10.60 5.34 8.62 一 一 7.18 
砂粒 25.62 10.50 6.12 5.50 1.37 3.40 一 一 30.72 
壤 粒 28.72 13.43 5.14 5.40 0.84 3.00 一 一 24.65 
黏 粒 31.75 20.44 8.46 6.94 3.43 2.62 一 一 12.01 
光谱 主 份 1 21.40 14.41 8 87 13.73 1.96 10.22 一 一 13.19 
光谱 主 份 2 24.74 19.07 10.21 10.01 6.13 2.03 一 一 14.92 
样 带 2 有 机 质 30.61 872 8.57 9.41 16.91 0.58 0.48 0.22 8.94 
Transect 2 高 程 1.47 1.39 1.20 j:71 2.33 1.35 1.16 0.63 72.83 
一 坡度 17.42 9.36 10.04 9.48 13.26 3.23 2.08 2.59 5.20 
es 湿度 指数 25.55 15.69 12.67 9.74 5.80 3.41 0.77 0.02 8.68 
2 容重 37387 11.58 9.31 10.49 3.76 0.39 0.63 0.23 9.49 
砂粒 19.99 7.69 7.36 2.99 19.32 2.16 0.46 1.33 20.36 
款 粒 22.17 13.94 8.05 4.17 21.16 0.71 0.12 0.47 10.40 
黏 粒 30.71 12.92 9.45 5.70 5.91 3.39 0.76 1.41 15.83 
光谱 主 份 29.29 14.60 13.61 12.80 9.44 1.46 0.30 0.12 5.68 
光谱 主 份 2 48.67 13.92 8.32 8.73 8.73 0.58 0.04 0.32 1.41 
样 带 3 有 机 质 17.62 6.14 3.35 5.35 12.53 20.46 4.98 一 2.44 
Transect 3 高 程 0.54 0.87 1.31 0.87 2.37 21.07 2.37 一 70.04 
坡度 30.20 9.87 11.58 5.58 5.76 11.34 0.86 一 6.20 
湿度 指数 39.22 13.95 8.57 6.66 4.27 5.02 0.08 一 2.21 
容重 26.22 12.71 8.38 10.66 9.32 1.73 0.83 一 8.03 
砂粒 25.96 13.41 10.12 7.82 7.67 2.17 1.55 一 20.08 
壤 粒 25.89 15.90 12.64 5.50 6.93 2.68 0.59 一 11.32 
黏 粒 34.23 10.91 14.70 6.46 6.76 1.77 1.20 一 7.81 
光谱 主 份 1 16.42 5.13 4.61 5.98 19.57 16.47 3:27 一 20.02 
光谱 主 份 2 21.98 7.03 7.17 6.90 11.34 20.45 2.38 一 9.30 


通过 观 


察 


各 影响 因子 的 IMF 方差 百分比 发 现 ,高 程 除 在 第 三 条 样 带 IMF6 处 值 (21.07% ) 较 大 外 ,其 余 主 
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要 分 布 在 残 差 部 分 (3 条 样 带 分 别 为 75.93% ,72.83% 和 70.04% ) ,表明 高 程 的 主要 表征 尺度 比 该 实验 提取 的 
尺度 大 。 在 样 带 1 中 , 除 高 程 外 的 其 他 影响 因子 的 表征 尺度 主要 在 IMF1 和 2 处 , 即 小 尺度 1011 和 1725 m 
处 。 样 带 2 中 ,坡度 .砂粒 和 壤 粒 含量 的 表征 尺度 主要 在 IMF1 和 5 处 (尺度 982 和 8573 m) ,其 余 影响 因子 主 
要 在 IMF1 和 2 处 (尺度 982 和 1530 m)。 样 带 3 中 , 除 坡度 (IMF1.3 和 6 处 ) re 和 3 处 ) .光谱 主 份 
1 和 2(IMF1.5 和 6 处 ) 外 ,其 余 影响 因子 主要 在 IMF1 和 2 处 (960 和 1504 m)。 结果 表明 , 除 高 程 外 相关 
因子 的 变异 性 主要 发 生 在 小 尺度 处 。 

与 采样 尺度 上 的 相关 性 不 同 ,土壤 有 机 质 与 影响 因子 多 尺度 相关 性 见 表 5$。 除 样 带 3 中 IMF2 外 , 商 程 与 
土壤 有 机 质 的 相关 性 主要 表现 在 大 尺度 , 即 样 带 1 为 IMF5 .6 和 残 差 ; 样 带 2 为 IMF6 ,7.8 和 残 差 ; 样 带 3 为 
IMF6 7 和 残 差 。 这 些 结果 表明 ,盆地 内 高 程 对 土壤 有 机 质 的 影响 主要 表现 在 大 尺度 约 10000 12000 .23000 m 
处 。 坡 度 与 有 机 质 关 系 在 样 带 1 中 较 弱 ,主要 表现 在 大 尺度 IMF6 和 残 差 处 ,在 样 带 3 中 主要 表现 在 JIMF2 .5、 
6、7 和 残 差 处 ;而 与 样 带 2 中 有 机 质 关 系 最 为 显著 , 即 除 IMF1 和 3 外 ,其 余 均 显著 相关 4 地 形 湿度 指数 与 土 
培 有 机 质 的 空间 多 尺度 关系 表明 ,在 盆地 中 部 两 者 的 关系 最 明显 , 主要 表现 在 尺度 98242588 4487 .8573 、 
10901 11591 m 处 。 盆 地 下 部 两 者 的 关系 较 明 显 ,这 与 采样 太 度 上 两 者 的 关系 不 一 致 。 羽 壤 容 重 与 有 机 质 的 
多 尺度 关系 表明 ,在 盆地 上 .中 部 小 尺度 约 1000 和 1500 m 处 ,中 下 部 尺度 约 3000 mR 下 部 尺度 约 


表 5 基于 MEMD 的 不 同 IMF 和 残 差 的 土壤 有 机 质 与 影响 因 巴 的 相关 系数 
Table S ”Correlation coefficients between SOM and influencing factors for each IMF and residue based on MEMD 


样 带 本 征 模 高 程 坡度 这 度 指 数 。 容量 动 奏 和 站 和 Es 
Transedt 隐 数 joation Slope ee Bu Sd Silt Clay 主 份 1 主 份 2 
IMF index density PC1 PC2 
样 带 1 IMF1 一 0.06 0.14 0.18 —0.308™ -0.20 0.24 0.01 0.58 —0.14 
Transect 1 IMF2 —0.08 —0.02 0.03 -0.4087™ 二 0.33 0.29 一 0.12 0.61 一 0.20” 
IMF3 -0.16 0.02 0.41 0.12 -0.06 0.37 一 0.27 0.32 0.08 
IMF4 0.06 -0.13 -0.34”” 0-0.20™ 0.30”” -0.10 一 0.36 一 0.90 一 0.71 
IMFS 0.35 0.05 0.00 0.54 一 0.08 0.32 一 -0.11 0.43 -0.10 
IMF6 -0.83*” -0.207 一 0.17 一 0.77 -0.09 0.40 一 0.353 0.77 一 一 0.07 
残 差 -0.41”* -0.27 全 0.39 一 0.57 一 0.26”” -0.21” 一 0.34 0.86 一 0.33 
样 带 2 IMF1 一 0.04 0.12 -0.24”” -0.33”” —0.34™”™ 0.28™™ 0.19” 0.51 一 0.09 
Transect 2 IMF2 0.06 0.41 0.02 =0.32** 一 0.38 一 0.49 一 —0.14 0.33 0.05 
IMF3 0;12 0.09 -0.21”™ 一 0.34 一 0.41 0.33 0.23 0.58 一 一 0.39 一 
IMRF4 0.05 0.26™* 0.57™™ 0.03 一 0.31 0.40 一 0.05 0.72 0.07 
IMFS =0.03 0.68™* 0.38 一 0.10 一 0.78 0.85 0.43 0.58 一 一 0.60 一 
IMF6 -0.25”* -0.41** 0.30”* -0.62™* —0.04 0.48™™ —0.23** 0.73™* 0.46™* 
TMF7 -0.90”” -0.74™* 一 0.54” -0.95™” 0.03 0.50 一 0.30 0.88 一 0.95 一 
IMF8 一 0.49 0.41 0.10 一 1.00 0.64”*” -0.57 一 一 0.69 一 0.96 一 一 0.46 一 
己 差 -076** 091% -036% -0.99* -098 096 099… 0.98" -0.41™ 
样 带 3 IMF1 一 0.07 0.09 0.03 -0.12 -0.29™" 0.38 一 0.03 0.53 0.06 
‘Transect 3 IMF2 -0.20° 一 0.23 一 一 0.07 0.08 一 0.353 0.40 一 0.00 0.44 一 -0.13 
IMF3 —0.11 -0.06 0.19™ 一 0.43 一 0.27 一 0.30 -0.01 0.57 0.28 一 
IMF4 0.04 0.04 0.51*” -0.35”* -0.26™” 0.19” 0.22 和 0.62 一 0.26 一 
IMFS 0.00 0.29™* —0.15 一 0.67 一 0.24 人 227 天 0.06 0.81 一 0.33 
IMF6 -0.75*” -0.61™* -0.50”” -0.83”* 0.11 0.37 一 一 0.73 0.80 一 0.21“ 
IMF7 一 0.46 ”一 0.99 一 0.88” -0.25 一 一 0.31 0.81 0.17 0.90 一 一 0.84 一 
残 差 0.40”” -0.43*” -0.73 -0.19” 0.93 ” 一 0.90 一 一 0.93 一 0.82 一 一 0.97 一 


* 显著 性 水 平 为 P<0.05; * * 显著 性 水 平 为 P <0.01 
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4500 和 6700 m 处 ,上 部 尺度 约 5400 和 9500 m 处 及 在 整个 盆地 内 大 尺度 > 10000 m 处 ,两 者 显 车 相关。 土壤 
质地 (包括 砂粒 、 壤 粒 和 黏 粒 ) 与 有 机 质 的 多 尺度 关系 表明 , 塘 粒 含量 对 土壤 有 机 质 的 多 尺度 关系 最 显著 , 即 
除 样 带 1 中 IMF4 外 ,二 者 在 所 有 表征 尺度 均 显 著 相 关 。 光 谱 主 份 与 有 机 质 的 多 尺度 关系 表明 ,在 研究 样 带 中 
的 所 有 表征 尺度 上 ,光谱 主 份 1 均 与 有 机 质 含量 显著 相关 ,而 光谱 主 份 2 的 相关 性 明显 弱 于 光谱 主 份 1。 总 
之 ,相关 因子 对 土壤 有 机 质 的 影响 随 尺度 增 大 而 呈 增 大 趋势 。 另 外 ,对 于 残 差 部 分 ,土壤 有 机 质 与 影响 因子 的 
相关 性 均 达 到 0.05 的 显著 性 水 平 ,表明 被 选 定 的 影响 因子 与 土壤 有 机 质 存 在 一 定 关系 。 
2.3 ”基于 与 影响 因子 多 尺度 关系 的 土壤 有 机 质 预测 

如 表 6 所 示 ,本 研究 采用 逐步 多 元 回归 获取 了 每 个 IMF 土壤 有 机 质 的 预测 模型 。 从 中 可 以 看 出 , 除 样 带 
3 中 IMF2 外 ,所 有 预测 模型 的 可 调整 R 均 在 0.44 一 1.00 之 间 变 化 , 且 随 IMF 增 大 而 增 大 。 可 调整 R 的 变化 
趋势 表明 ,尺度 越 大 ,对 土壤 有 机 质 的 预测 精度 越 高 , 且 模型 中 所 选 因 子 对 有 机 质 的 影响 更 强烈 。 攻 于 所 有 
IMF 和 残 差 的 有 机 质 预 测 结果 ,本 人 研究 采用 逐步 多 元 回归 法 进行 采样 尺度 上 土壤 有 机 质 的 预测 ;结果 列 于 表 
7。 结 果 表 明 : 采 样 尺 度 上 各 样 带 中 土壤 有 机 质 预测 值 和 实测 值 的 尼 分 别 为 0.90.0.86 和 和 0557 显著 高 于 采样 
尺度 上 直接 利用 逐步 多 元 回归 拟 合 的 结果 (0.52 .0.56 和 0.54) 。 通 过 MEMD 方法 获取 的 有 RMSE 和 NRMSE 显 
著 低 于 直接 逐步 多 元 回归 的 预测 结果 ,而 RPD 显著 高 于 直接 逐步 多 元 回归 的 预测 结果 。 基 于 上 述评 价 参数 
的 比较 ,可 以 得 出 采用 MEMD 对 土壤 有 机 质 的 预测 精度 高 于 基于 原始 数据 的 简单 逐步 多 元 回归 预测 的 结论 。 


表 6 基于 经 验 模 态 分 解 的 每 个 IMF 和 残 差 的 土壤 有 机 质 逐 步 多 元 回归 预测 异型 和 回归 统计 特征 ( 玉 值 和 可 调整 R?) 
Table 6 ” Predictive equations and regression statistics ( F-value and adjusted R?) fomSOM for each IMF and residue using stepwise multiple 
linear regression based on multivariate empirical mode decomposition 


样 带 公式 


IMF F 2 
Transects Function A 
样 带 1 IMF 1 0.03+21.31(0.43)Silt+3.77(0.68)PC1-4.94(082) PC2 42.0 0.52 


0.01+1.06( 0.17) TWI-12.01(0.19) BD-27.28( 0.46) Sand—9.66( 0.26) Clay+ 
3.63(0.53)PC1-21.01(0.42) PC2 


0.05+0.07( 0.16) Elevation#OW26( 0W13) Slope+ 1.15(0.24)TWI +9.74( 0.44) BD+ 
22.34( 0.37) Sit-18.42(031) Clay+3.63( 0.81)PCI1 


IMF 4 —0.06+0.03( 0.09) Eevation—0.34( 0.17) Slope+4.96( 0.07) Silt+0.69( 0.90) PC1 142.6 0.83 
—0.03+0.23( 0,48)Elevationt1.46( 0.70) Slope+1.55(0.27)TWI-7.08( 0.48) BD+ 


Transect 1 IMF 2 


IMF 3 53:2 0.77 


IMF 5 1808. 0. 
58.61(0.59)Silt+6.17(0%10) Clay+4.54(0.72)PC1-0.21(0.24)PC2 © 2 
IMF 6 0.01+0.46(0s98) Elevation+3.62(0.75)Slope+6.44(0.18)TWI-3.75(0.25)BD- Gd 1.00 
157.71(2.27) Sand—204.82( 1.49)Clay+7.41( 1.38)PC1+ 7.49(0.07)PC2 
区 -1007.72+I161(2.02)Slope+138.15(1.38)BD+19.51(2.44)PC1+ 
避 394573. 1.00 
残 关 251.95(4.02)PC2 8 
样 带 2 IMF@ 0.07+0.49(0.14) Slope—1.41(0.27)TWI-10.56( 0.26) BD+6.77(0.12)Silt+ 有 
Transect 2 332(0.51)PC1 
Mr 0.03—0.16(0.11) Elevation—1.03(0.41) Slope—7.91( 0.20) BD-3.33(0.37) Sand- 2 0.59 
29.89(0.59)Clay+1.95(0.4)PC1+4.66(0.15)PC2 
.08—0. svation—1. = .13 二 1S， 
MIS 0.08—0.58(0.30) Elevation-1.12(0.29)TWI-S.67(0.13)BD-15.07(0.26 ) Clay+4.62 48.8 说 


(0.92)PC1-17.81(0.44)PC2 


-0.07+0.71(0.28) Slope+1.71(0.36)TWI-2.56(0.06)BD+1.42(0.16) Silt+ ee 0 
4.30(0.79)PC1-13.93(0.34)PC2 


-0.01-0.55(0.51) Elevation-0.52(0.18) Slope+2.67( 0.33) TWI+30.02( 0.72) Silt- 


IMF 5 418.7 0.95 
28.44(0.51) PC2 2 
0.33(0.85) Elevation+0.57( 0.52) Slope-2.42(1.22)TWI-37.26(0.68) BD +18.02 

6 (0.71)Clay+9.25(2.32)PC1-16.24(0.41)PC2 10392.0 1.00 
0.16( 0.30) Elevation-2.73(0.72)TWI-11.86(0.34)Sand-16.14(0.17) Sit+11.10 
IMF 7 12571476.7 1.00 
(1.39)PC1-33.83(0.25)PC2 8 
IMF 8 -0.03(0.18) Elevation+1.50(0.10)TWI-46.05(1.10)BD+11.81(0.36)PC2 1525552.4 1.00 
残 差 —32.33+1.41( 0.42) Slope+32.04( 0.69) Clay+17.29(0.17)PC2 816066.5 1.00 
样 带 3 IMF 1 -0.06-23.23(0.43)Sand—19.10(0.27) Clay+3.81(0.61)PC1-11.97(0.27)PC2 25.4 0.44 
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续 表 
涯 带 八 式 
样 带 让 公式 本 Rp2 
Transects Function 
Transect 3 IMF 2 —0.52(0.16)TWI+13.43(0.27)Silt+3.26(0.50) PC1-14.41(0.31)PC2 19.5 0.38 
IMF 3 0.0150.11 (027) levationt 0:44( 0.18 Slgpet 03 O33 WI BD i 0.49 
11.81(0.28)Silt+2.54(0.50) PC1-6.04(0.18)PC2 
0.20+2.98(0.66)TWI+12.61(0.26) BD+4.53(0.31)Silt-17.50(0.19) Clay+ 
.3 
IMF 4 5.77(1.01) PC1 216 0.89 
IMF 5 0.05—0.16(0.45) Elevation+2.91(0.44) Slope+1.01(0.12)TWI-15.34(0.19) BD- 323.9 (6 95 
8.10(0.07) Silt-47.86(0.35 ) Clay+2.52( 0.53)PC1-32.90( 0.64)PC2 
-0.4-0.24(2.08)Elevation+6.90(1.14) Slope-13.72(1.36)TWI-18.46(0.08)BD- 
IMF 2766. 0. 
0 48.27( 0.22) Silt+243.41(0.72 ) Clay > 2 
: —5—0. svation—2.31(0. 一] —78. 
IMF 7 7.91E-5-0.03(0.15) Elevation—2.31(0.21)Slope-1.16(0.03)TWI-78.23( 0.46) 148763491 1.00 


BD+3.64(0.49) PC1-16.36( 0.23) PC2 
残 差 117.97-0.73( 0.26) Slope+46.93( 0.84) Silt-43.50( 1.72) PC2 1452867 1.00 
”levation: 高 程 ;Slope: 坡 度 ;TWI; 地 形 湿度 指数 ,topographic wetness index; BD:; 容重 ,bulk density; Sand ;砂粒 ; Siib 坟 入 ;Clay: 条 粒 ;PC1; 光谱 
主 份 1; PC2 :光谱 主 份 2; 括号 中 的 数据 表示 标准 偏 回 归 系 数 


表 7 评价 土壤 有 机 质 预测 精度 的 相关 参数 
Table 7 Values of evaluation parameters used to assess the overallSOM prediction accuracy 


预测 方法 性 开 


Predicting methods Transects ue > a R 
逐步 多 元 回归 样 带 1 5.34 0433 1.45 0.52 
Stepwise multiple regression 样 带 2 3.81 0.22 1.51 0.56 
样 带 3 5.45 0.27 1.48 0.54 
经 验 模 态 分 解 样 带 3.35 0.21 2:32 0.90 
Multivariate empirical mode decomposition 样 带 2.90 0.17 1.98 0.86 
样 带 3 3.36 0.17 2.40 0.91 


RMSE : 均 方 根 误差 ,root mean squared error; NRMSE: 标准 化 均 方 根 误差 , normalized root mean square error; RPD : 相对 预测 偏差 ,residual 


prediction deviation 


采样 尺度 上 土壤 有 机 质 实际 值 与 单个 TMF 或 残 差 处 预测 值 及 变量 预测 有 机 质 值 的 相关 系数 如 图 3 所 
示 。 该 相关 系数 的 比较 可 反映 每 和 IMF 相对 于 采样 尺度 上 有 机 质 预 测 结果 的 重要 性 程度 。 结 果 表 明 , 样 带 1 
中 有 机 质 预测 的 主要 贡献 者 是 TMF6; 样 带 2 中 主要 贡献 者 是 IMF5 ; 样 带 3 中 主要 贡献 者 是 IMF6。 即 盆地 上 
部 尺度 12375 m 处 中 部 尺度 8573 m 处 和 下 部 尺度 11806 m 处 对 采样 尺度 上 有 机 质 的 预测 起 主要 作用 。 另 
外 ,可 能 代表 更 大 尺度 的 残 差 在 样 带 1 中 的 比重 较 大 ,其 次 是 样 带 2 ,最 后 是 样 华 3 ,表明 被 分 解 的 IMF 能 够 很 
好 解释 土壤 有 机 质变 异性 的 顺序 为 盆地 的 下 部 > 中 部 > 上 部 。 

土壤 和 有机质 的 变异 性 主要 是 土壤 本 身 异 质 性 地 形 植被 人 类 活动 等 综合 因素 影响 的 结果 。 不 同 尺度 的 
某 一 因 予 对 有 枯 质 预测 值 与 采样 尺度 上 实测 值 的 相关 系数 如 图 3 所 示 。 该 结果 可 表示 每 一 因子 对 采样 尺度 
上 土壤 有 机 质 预 测 的 相对 重要 性 。3 条 样 带 中 光谱 主 份 1 (综合 因子 ) 是 有 机 质 预 测 中 最 重要 的 影响 因子 。 
此 四 Si 样 带 1 中 容重 、 样 带 2 中 砂粒 和 壤 粒 含量 样 带 3 中 地 形 湿 度 指数 .坡度 和 容重 对 其 有 机 质 预 测 也 起 了 
重要 作用 。 总 之 ,地 形 因 子 .土壤 容重 和 光谱 主 份 对 样 带 3 中 有 机 质 的 影响 最 为 明显 ,土壤 质地 对 样 带 2 中 有 
机 质 的 影响 最 为 明显 , 且 影 响 顺 序 为 中 部 > 下 部 > 上 部 ,而 样 带 1 中 土壤 有 机 质 受 地 形 因 子 土壤 容重 和 质地 
的 影响 程度 最 弱 。 

与 传统 回归 分 析 相 比 ( 表 2) ,MEMD 可 捕捉 到 样 带 1 中 容重 、 样 带 2 中 砂粒 及 样 带 3 中 地 形 湿 度 指数 和 
坡度 对 有 机 质 的 影响 。 同 时 ,回归 方程 中 标准 回归 系数 也 可 表示 这 些 影响 因子 在 每 一 尺度 上 对 有 机 质 预 测 的 
相对 重要 性 ( 表 6) 。 在 样 带 1 中 ,光谱 主 份 1 对 所 有 表征 尺度 上 土壤 有 机 质 预测 起 显著 影响 。 土 壤 质地 在 小 
尺度 (大 3078 m) 和 大 尺度 (三 9535 上) 处 对 有 机 质 预测 起 显著 影响 。 同 时 ,地 形 因子 对 有 机 质 预测 的 影响 随 
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影响 因子 Influeneing factors 


图 3 采样 尺度 上 土壤 有 机 质 实际 值 与 单个 IMF 或 残 差 处 预测 值 及 变量 预测 有 机 质 值 的 相关 系数 
Fig.3 “Coefficient of determination between SOM at the measurement scale and predicted IMFs (residue) or total SOM predicted by 


each variable 


尺度 增 大 而 增强 ,如 坡度 因子 在 IMF3、4 和 人 @ 处 系数 分 为 0.13 .0.17 .0.70 和 0.75。 土 壤 容重 对 IMF2 .3 .5 .6 
和 残 差 处 有 机 质 预 测 都 有 显著 影响 。4 因 此》 容重 对 样 带 1 采样 尺度 上 有 机 质 预测 的 贡献 主要 体现 在 这 几 个 尺 
度 上 的 综合 作用 。 

在 样 带 2 中 ,由 于 IMF7 和 8 中 有 机 质 预测 值 与 采样 尺度 上 实测 值 的 相关 系数 不 显著 , 故 可 忽略 (图 3)。 
从 IMF1 到 6 处 ,地 形 因子 对 有 有 有机质 预测 的 影响 逐渐 增强 。 土 壤 容 重 和 质地 对 不 同 尺度 有 机 质 预测 的 影响 无 
明显 规律 。 砂 粒 含 量 对 采样 尺度 上 有 机 质 预测 的 影响 主要 体现 在 IMF2 处 。 在 样 带 3 中 ,地形 因子 对 有 机 质 
预测 的 影响 随 尺 度 的 增 大 而 增强 ,上 且 在 IMF6 处 达到 最 大 值 ,然后 随 尺 度 增 大 而 减 小 , 且 湿度 指数 对 样 带 3 中 
有 机 质 预 测 的 影响 最 为 明显 。 土 壤 质地 在 其 小 尺度 IMF1 和 大 尺度 IMF6 处 影响 最 强 。 容 重 在 IMF3 .4 和 7 
处 对 有 机 质 痊 测 寡 稳 是 影响 ,其 余 尺度 不 明显 。 综 上 ,与 传统 回归 分 析 相 比 , MEMD 方法 在 某 些 表征 尺度 可 
捕捉 到 相关 因 巴 对 有 机 质 分 布 的 影响 , 故 其 预测 精度 高 于 采用 传统 回归 分 析 方 法 。 

韦 庆 有 机 质 实测 值 和 预测 值 的 局 部 小 波 如 图 4 所 示 。 与 MEMD 相 比 ,逐步 多 元 回归 预测 误差 在 样 带 1 
中 主要 是 在 位 置 19.8 一 33.0 km 处 局 部 高 佑 ,造成 尺度 1 一 2 和 4 km 处 局 部 方差 的 增 大 ; 样 带 2 中 ,主要 是 在 
位 置 49.8 一 33.0 km 处 局 部 低估 ,造成 尺度 4 一 8 km 处 局 部 方差 的 减 小 ; 样 带 3 中 ,主要 是 在 位 置 19.8 一 36.3 
km 处 局 部 高 估 ,造成 8 km 处 局 部 方差 的 增 大 。 


本 研究 采用 MEMD 法 将 太原 盆地 不 同位 置 的 土壤 有 机 质 空 间 序列 分 解 为 不 同 的 表征 尺度 ,发现 主 要 表 
征 尺 度 在 1000 m 处 。 因 此 ,在 利用 栅 格 数据 存储 该 盆地 内 土壤 属性 时 ,建议 将 表层 土壤 有 机 质 最 佳 格 网 宽度 
设置 为 1000 m, 以 便 捕捉 到 有 机 质 较 大 的 空间 变异 性 。 伟 地 上 部 残 差 的 方差 百分比 较 大 (24.19%) ,可 能 是 需 
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图 4 土壤 有 机 质 实测 值 和 逐步 多 元 回归 、MEMD 法 预测 值 的 局 部 小 波 图 
Fig.4 Local wavelet spectrum of measured and predicted SOM by SMLR or MEMD 


设计 更 长 的 采样 带 才 可 利用 MEMD 识别 的 更 大 尺度 。 同 时 , 侈 地 上 部 的 表征 尺度 与 中 、 下 部 差异 较 大 ,可 能 
是 由 于 盆地 上 部 距 太原 和 晋 中 市 等 大 .中 型 城厢 较 近 ,人 为 活动 造成 十 地 利用 破碎 。 此 外 ,有 机 质 空间 序列 在 
垂直 河流 方向 上 的 主要 表征 尺度 沿 汾 河流 域 放 向 表现 呈 分 散 状 态 ,表明 沿 河流 方向 的 相关 因子 对 有 机 质 空间 
分 布 影响 的 主导 性 增强 。 

在 采样 太 度 和 空间 多 尺度 关系 -盆地 中 部 土壤 有 机 质 分 布 与 地 形 因 子 容重 和 土壤 质地 的 相关 性 最 强 ， 
这 可 能 是 由 该 区 处 于 侵蚀 平原 (_B 部 ) 与 陷落 盆地 (下 部 ) 的 交界 处 的 特殊 地 质 条 件 造 成 。 在 采样 尺度 上 ， 
盆地 中 部 有 机 质 与 光谱 主 份 贡 的 相关 性 最 弱 ,可 能 是 由 于 中 部 区 域 土壤 有 机 质 的 变异 性 较 小 ,导致 有 机 质 与 
土壤 光谱 之 间 的 相关 性 减弱 "" 。 另 外 ,有 机 质 与 某 些 因子 在 采样 尺度 不 存在 显著 相关 ,但 是 二 者 的 空间 多 尺 
度 关 系 在 某 些 表征 尺度 呈 显 著 相 关 , 表 明 土 壤 有 机 质 与 影响 因子 的 空间 多 尺度 关系 能 获取 更 多 信息 。 同 时 ， 
rn A s 间 多 尺度 关系 的 有 机 质 预测 精度 显著 高 于 采样 尺度 上 直接 利用 逐步 多 元 回归 分 析 

结果 , 进 合 步 表明 单一 尺度 (采样 尺度 ) 的 简单 分 析 难 以 揭示 两 者 的 复杂 关系 。 

MEMD 法 通用 于 分 析 非 线性 . 非 平 稳 数 据 序列 的 空间 多 尺度 关系 "的 。 由 于 该 理论 建立 时 间 尚 短 ,在 土壤 
属性 的 空间 多 尺度 研究 中 还 未 广泛 应 用 。 同 时 ,MEMD 法 也 存在 一 定 缺 陷 。Hu 和 Si 研究 结果 表明 ,每 条 样 
带 的 所 有 IMF 和 残 差 的 方差 百分比 之 和 并 不 等 于 100%!231 。 此 外 , 与 小 波 分 析 不 同 ,MEMD 法 会 针对 不 同 采 
样 带 分 解 出 不 同 的 表征 尺度 ,特定 样 带 的 研究 结果 不 能 推广 。 因 此 ,MEMD 法 重 在 分 析 其 发 生 过 程 ,在 结 
验证 方面 存在 缺陷 。 


4 结论 


本 文采 用 MEMD 法 分 析 了 太原 盆地 土壤 有 机 质 与 影响 因子 的 空间 多 尺度 关系 ,主要 结论 如 下 : 
(1) 盆 地 上 部 主要 表征 尺度 为 1011 和 1725 m ,中 部 为 982 和 8573 m ,下 部 为 960 .6753 和 11806 m。 整 
个 盆地 内 ,尺度 约 1000 m 处 是 土壤 有 机 质 的 主要 表征 尺度 , 且 盆 地 内 垂直 河流 方向 的 有 机 质 序 列 主要 表征 尺 
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度 沿 河流 方向 表现 分 散 。 
(2) 土壤 有 机 质 与 影响 因子 在 采样 尺度 和 MEMD 空间 多 斥 度 的 相关 性 顺序 为 :盆地 中 部 > 下 部 > 上 部 。 
(3) 在 3 种 景观 样 带 上 ,光谱 主 份 1 与 有 机 质 的 相关 性 均 显 著 。 其 次 ,盆地 上 部 的 容重 .中 部 的 砂粒 和 下 
部 的 地 形 湿 度 指 数 对 其 影响 较 明 显 ,而 在 采样 尺度 上 念 地 下 部 二 者 的 关系 并 不 显著 。 因 此 ,单一 尺度 分 析 不 
能 够 全 面 揭示 土壤 有 机 质 与 相关 因子 在 所 有 空间 扩 度 上 的 复杂 关系 ,而 MEMD 法 对 有 机 质 的 预测 精度 要 显 
著 高 于 直接 利用 逐步 多 元 回归 分 析 。 


致谢 :感谢 郭 颖 . 王 鹏 和 徐 振 对 样品 采集 和 测试 提供 的 帮助 。 
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